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Le domaine d'application de Tinvention est celui de la prospectlon 
sismique. LMnvention concerne plus particulierement le traitement des 
traces sismiques d'une collection en point milieu commun. 

Linvention concerne plus pr6cisement un precede de determination 
5 des parametres de Vitesse V et d'anellipticite rj necessaires pour realiser 
des traitements comprenant une correction d'obliquit^ des traces 
sismiques. 

La prospectlon sismique consiste d'une maniere generate. §i §mettre 
dans le sous-sol, k Taide d'une ou plusieurs sources sismiques, des ondes 

10 sismiques, a enregistrer en surface, en fonction du temps, des donn^es 
sismiques correspondant aux ondes sismiques r§flechies sur les interfaces 
g^ologiques du sous-sol (encore appel6es' reflecteurs) a Taide de 
r§cepteurs (encore appeI6s g^ophones ou hydrophones selon que Ton 
prospecte a terre ou en mer) puis a traiter ces donnees pour en extraire des 

15 informations utiles quant k la geologie du sous-soL 

On appelle trace sismique Tenregistrement de T^nergie sismique 
realise par chaque recepteur durant I'acquisition des donnees. 

Une technique classique de prospectlon sismique est la couverture 
multiple pour laquelie sources et rdcepteurs sont agences de telle sorte 

20 qu'un m§me point milieu (c'est-gi-dire le point a egale distance entre la 
source et le recepteur d'une trace considSree) regroupe plusieurs traces 
sismiques. 

Si les traces sismiques contiennent des informations utiles quant aux 
reflexions sismiques et S la geologie du sous-sol, elles contiennent 
25 egalement des composantes de bruit 

L'un des premiers objectifs du traitement des donnees sismiques est 
d'eliminer, ou tout du moins d'attenuer, ces composantes non desirees de 
bruit de telle sorte que les informations utiles puissent etre clairement 
identifi6es et interpr6tees. 
30 Une methode classiquement utilis6e afin d'attenuer ces 

composantes de bruit est la collection en point milieu commun (ou 
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collection CMP selon I'acronyme de I'expression anglo-saxonne Common 
Midpoint). Les traces disposant d'un m§me point milieu sont alore 
regroup6es selon la distance s§parant source et r6cepteur (appelee deport 
ou « offset » selon la temninologle anglo-saxonne). 
5 D'une manidre g6n6rale. la repr6sentation en image des donnees 

sismlques n§cessite la mise en oeuvre d'un traitement comprenant : 

• une operation dite TZO (selon I'acronyme de rexpression anglo- 
saxonne Transform to Zero Offset) visant a compenser I'effet 
d'obliquit6 des trajets en ramenant les temps d'amV§e des reflexions S 
10 ceux de traces a deport nul, 

. ainsi qu'une operation de migration visant ^ restituer les formes 
correctes des interfaces g§ologiques. 
Si ces operations de TZO et de migration sont g6n6ralement 
realisees successivement. elles peuvent aussi §tre realis§es conjointement. 
15 Cela est en particulier le cas lorsque est realis6e une migration temps avant 
sommation (migration PSTM selon I'acronyme de I'expression anglo- 
saxonne Pre-StackTime Migration). 

De manldre simplffiee. l'op6ration de TZO permet de simuler 
I'acqulsltfon des donnees sismlques par des sources et r§cepteurs disposes 
20 au point milieu commun. 

L'objectif est d'additionner les enregistrements illuminant le m§me . 
point du sous-sol afirvd'augmenter le rapport signal sur bmit et le rapport 
reflexions primaire sur reflexions secondaires, et de beneficier ainsi des 
atouts de la « couverture multiple ». 

Afin de fabriquer I'image § d§port nul, une m6thode dite de con-ection 
d'obliquite ou NMO (correspondant ^ I'acronyme de I'expression anglo- 
saxonne Normal Move Out) est mise en oeuvre. 

Si on fait I'hypothdse d'un sous-sol stratlfi§ horizontalement et sans 
variation Iat6rale des vitesses de propagation, on montre que les 
enregistrements ayant la propriet§ d'eclairer le m§me point du sous-sol sont 
ceux ayant le m§me point milieu. 
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Cependant I'image d'une reflexion dans le sous-sol arrive & des 
temps variables selon le deport. Afin d'additionner les reflexions, il faut done 
prealablement corriger les differents enregistrements pour les ramener tous 
d une r6ference commune, celle de deport nul, 
5 Historlquement, la correction d'obliquit§ repose sur un module 

particuli^rement simple du sous-sol : le moddle homogdne avec reflecteurs 
horizontaux. 

Dans ce module, les reflexions associ^s S chacun des rSflecteurs 
du sous-sol, s'alignent th6oriquement le long d'hyperboles, encore 
10 appel6es Indicatrices, centrSes a la vertlcale du point milieu. 

Le temps d'arriv§e d'une reflexion est alors une fonction 
hyperbollque du deport source-recepteur, le temps le plus court §tant celui 
obtenu & deport nul. 

De maniere & r§aliser la sommation des enregistrements de ohaque 
15 collection, la correction NMO redresse les hyperboles afin de les amener '4 
theorlquement ^ I'horizontal. 

La correction NMO est alors r§alis§e en se basant sur I'^quation > 
hyperbollque suivante du temps t de propagation apr§s reflexion, associe a 
un couple source - r6cepteur de deport x : 

dans laquelle to repr6sente le temps de propagation pour un deport nul et V 
designe la Vitesse m(^enne de propagation des ondes dans le sous-sol. 

Le modele simpliste evoque ci-dessus s'appuie notamment sur des 
hypotheses telles que de faibles angles d'incidence et un milieu isotrope. 

25 Mais les hypotheses du modele simpliste sent trop rest'rictives pour 

decrire un milieu complexe, et ne peuvent en particulier s'appliquer a la 
propagation des ondes sismiques dans un milieu anisotrope (milieu dans 
lequel la Vitesse des ondes peut varier selon la direction de propagation). 
L'utilisation d'un modele moins simpliste est rendue necessaire, en 

30 particulier du fait de : 
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- rutilisation de longues flQtes pour racquisition en off-shore profond, ce 

qui conduit en particulier d acquerir des traces d large deport ; 
^ Tobservation d'anisotropfe dans les sediments de type arglieux. 

L'hypothese g6n6ralement acceptee conslste ^ mod6liser un milieu 
5 anisotrope comme un empilement de couches isotropes transversaiement 
disposant d'un axe de symetrie vertical On parte alors d'anisotropie VTI 
(acronyme de I'expression anglo-saxonne Vertical Transverse Isotropy), 

II a ainsi ete propose de determiner les corrections d'obliquite devant 
§tre realis6e en : 

10 - introduisant de rinhomogen6it§ verticale dans un modele de milieu 
VTI homogene, comme cela a 6t6 pr6sent6 dans le document 
«ALKHALIFAH. T. et TSVANKIN, I., 1995, Velocity analysis for 
transversely Isotropic media : Geophysics, 60, 1550-1566 » ; 
ou encore, en 

15 - introduisant de Tanlsotropie VTI dans un modele ^ couches 
isotropiques stratiflees, comme cela a 6te montr6 dans le document 
« SILIQI, R. et BOUSQUIE, N., 2000, Anelliptic time processing based 
on a shifted hyperbola approach, 70th Ann, Internal Mtg.: Soc. Of 
Expl.Geophys., 2245-2248 ». 
20 Cette seconde approche, combinant inhomogeneite verticale et 

anisotropie VTI pour foumir un nouveau modele du sous-sol, semble la 
meilleure dans la pluparl des cas reels etudies. 

L'equation suiv^nte de correction par hyperbole d^calee aneliiptique 
du temps t de propagation apr§s reflexion, associ^ ^ un couple source - 
25 r6cepteur de deport x, decoule de ce moddle : 

oCi V represente la Vitesse classiquement utilisee en sismique, 
correspondant ^ de faibles deports, et tj est un parametre, dit anellipticite. 
I! a egalement 6t6 propose dans le document « SUAUDEAU, E. et 
30 SILIQI, R., 2001, Anelliptic pre stack time migration, Annual international 
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Meeting, CSEG Expanded Abstracts » d'inclure la correction d'obliquite par 
liyperbole d^cal6e anelliptique dans {'Equation du temps de trajet utilisSe 
lors de ('operation de migration PSTiy4. 

Liquation de la migration PSTM s'exprime classiquement sous la 
forme d'une double hyperbole d^l^ anelliptique, somme de deux radnes 
carries (Equation DSQR, seion i'expression anglo-saxonne Double SQuare 
Root;. 

En tenant compte de raneliipticit§, I'expression de cette Equation en 
deport constant devient la suivante : 



10 t=- 



1+817 



^0 + 



i 



(1 + 877)^^ 



1 



1+8;? 



ou : 



- les parametres V et 77 sont ceux mentionn^s ci-dessus. 

- . Xm represente les coordonn§es des points milieux, 

~ x-Xm represente Touverture (ou « aperture ») de la migration, 
15 - h represente le demi-deport source - r6cepteur, 

- to represente le temps double a Touverture nulle de Top^rateur 

On notera que lorsque i'ouverture x-Xm de la migration est nulle, 
Tequation (1b) de correction PSTM devient Tequation (1a) de correction 
NMO. La correction d'obliquite NMO constitue ainsi un cas particulier de la 
20 migration PSTM : celu^e la migration PSTM d'ouverture nulle, 

Finalement, afin de reallser un traitement de donnSes sismique 
comprenant une correction d'obliquit6 qui tienne compte de rh§terogen§lte 
verticale et de Tanisotropie de type VTI, il est done necessaire de 
determiner les deux parametres de vitesse V et d'anellipticit6 77 

25 Uestimation desdits parametres V et 77 peut etre classiquement 

realisee en deux passes telles que : 



I 
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- au cours de ia premidre passe, la distribution des vitesses V le 
long de I'axe des temps est estim6e en n'utilisant que les donates 
^ deports proohes ; 

- au cours de la seconde passe, ranellipticit6 7 est estim6e. le long 
de I'axe des temps, en utilisant : 

o la distribution des vitesses d§termin6e au cours de la 

premiere passe, et 
o I'ensemble des donn6es (y compris celles § larges deports), 
li a 6galement 6t6 montre dans le document : 
« SILIQI. R., 2001, Technological leap In time processing focuses the data 
throughout anisotropic media: First Break, 19, n"1 1, 612-618 », 
qu'une estimation des paramdtres V et 7 peut §tre r6alls6e en une seule 
passe, au cours de laquelle on realise des analyses bispectrales permettant 
de pointer simultan§ment les deux parametres V et 7 le long de I'axe des 
15 temps, et cela en utilisant toutes les donn6es. 

Dependant, les analyses denses des parametres de correction 
d*obliquit§ sont pr6f§rablement r6alis6es lorsque I'^quation de correction ne 
depend plus du temps to (on parte alors de correction d'obliquite statique). 
Une con-ectipn statique pennet effectivement de d§caler, pour un 
20 deport donn6, I'ensemble des 6chantlllons constituant chacune des traces 
d'un m3me temps 6t 

. Ainsi, lorsqu'uige correction statique est r6alis6e, le nombre de 
calculs devant gtre effectues peut etre significativement reduit et le 
phenomene d'etirement des traces est ellmin6, ce qui rend vlables lesdites 
25 analyses denses. 

Seules les analyses de vitesses et d'anellipticit§ en deux passes 
peuvent aujourd'hui etre mises en oeuvre. en particulier gr^ce d des 
approximations paraboliques des reslduelles d'obliquit6. de manidre a 
obtenir un pointe dense desdits parametres V et 7. 
30 Dans ce cadre : 
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• on estime tout d'abord des vitesses r^siduelfes. en utilisant les donnees 
a faible deport, suite k une premiere estimation, des vitesses ; 
. • on estime ensuite I'anellipticit6, sur toutes les donn6es, en utilisant les 
mises d jour des vitesses r6alis6e pr^c^emment. 
5 Une loi d'effacement (ou « mute » selon la tenminologie angio- 

saxonne) doit ainsi §tre d6finie, pour I'estlmation des vitesses r6siduelles, 
afin de ne conserver parmi I'ensemble des donn6es que celles consider^es 
comme dtant §i faible deport. 

Or refficacite de I'analyse en deux passes est partlculi^rement 
10 sensible au choix d'une telle loi d'effacement. 

En outre, ce sont les donn§es ^ large deport qui sont principalement 
utilis§es pour restimation de ranelliptlcit§ rj . 

Mais, restimation de I'anellipticite r^alisSe pour les donndes d large 
deport n'est pas trds precise, si bien que ia correction r§alis6e est 
IS finalement Imprecise. 

Un but de I'lnvention est de pennettre de s'affranchir de ces 
limitations et inconv^nients. en proposant de realiser un traitement 
comprenant une correctton d'obilquite statique qui soit plus efflcace et plus 
precis. 

20 Plus pr^cisement, rinvention a pour objectif la detemiination dense 

des param^tres de Vitesse V et d'anellipticite tj en une seule passe utilisant 

I'ensemble des donn6es disponibles, c'est-a-dire en s'appuyant sur 

.•4 

I'ensemble de la gamme de deports. 

A cet effet, I'lnvention propose un precede de determination des 
25 paramdtres de vitesse V et d'anellipticite rj pour un traitement de traces 
sismlques d'une collection d point-milieu commun (CMP) comprenant une 
con-ection d'obilquite anelliptique, caract§rise en ce qu'il comporte : 

• une eiape preiiminaire de definition d'une pluralite de noeuds (dtn, 
r^), lesdits noeuds etant significatifs de parametres dtn et 
30 representant respectivement la correction d'obliquite' pour le deport 
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IS 



20 



25 



maximal et ie temps de propagation S d§port nul en coordonn6es 
hyperboliques. jadlte §tape pr6limrnalre etant suivie 

• pour chacun des noeuds (dtn, r^) d§flnis lors de l'6tape pr6liminaire. 
des etapes de correction d'obliquite des traces de ia collection CMP en 
fonction des valeurs desdits paramdtres dtn et au noeud consider^, et 
de calcul de la fonction de semblance associee h ladite correction 
d'obliquite pour Ie noeud consider^ ; et 

• pour cheque temps to de pointe. d'une etape comprenant ia 
determination du noeud (dtn(to).ro(to)) de semblance maximale. 

• et d'une 6tape finale de conversion des param^tres dtn(to) et rjto) 
de mani^re ^ obtenir les lols de vltesse V(to) et d'aneliipticite n{^). 

Selon un premier mode de realisation de I'invention, Ie traitement 
realise est une connection d'obliquite NMO statlque des traces sismiques. 

Selon un second mode de realisation de I'invention, Ie traitement 
realise est une migration PSTM des traces sismiques, ladite migration 
PSTM comprenant une con-ectlon d'obliquite PSTM statique desdites traces 
sismiques. 

Un aspect prefere. mais non limitatif, du precede selon I'invention 
conceme la definition des paramdtres dtn et relativement ^ la vltesse V 
et d ranellipticite /?. de manidre k assurer des con-ections d'obliquite 

statiques, selon dtn = + \{ 1 

D'autres aspects, buts et avantages de la presente invention 
apparaTtront mieux ^ la lecture de la description detailiee suivapte. faite en 
reference aux figures annexees sur lesquelies : 

- la figure 1a represente I'iiyperbole anelliptlque decaiee pemiettant de 
realiser la correction d'obliquite NMO et illustre la signification des 
parametres et dtn ; 
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- la figure 1b represente I'equation DSQR de la double hyperbole 
anelllptique d§calee de la migration PSTIVI et illustre la signification 
des paramdtres et dtn ; 

- la figure 2 illustre Veffet du parametre sur la courbe de r§flexion 
5 corrigee par la correction NMO ; 

r. les figures 3a et 3b repr6sentent le volume d'analyse (to. dtn, t„ ) dans 
lequel le point6 bispeotral selon I'invention des param^tres dtn et 
est rSalis^ ; 

- la figure 4 juxtapose chacune des approches (V, Van) et (dtn, t^) en 
10 representant leur panneau respectif d'analyse bispectrale ; 

- la figure 5 repr§sente la correspondance entre les deux paires (dtn, 
^0 ) et (V, Van) de paramdtres de correction d'obliquite ; 

- la figure 6 repr§sente une collection CMP de traces sismlques r§elles ' 
avant correction d'obliquit§ NMO ainsi que le point§ bispectral des 1' 

15 parametres (dtn, ) con-espondant d cette collection de traces ; ? 

- la figure 7 represente la fonction de semblance et les fonctlons de 
Vitesse et d'anelliptlcit6 d§duites du pointe des param§tres dtn et 
de la figure 6. 

- la figure 8a est un organigramme representant les 6tapes d'un premier 
20 mode de realisation particulier de I'invention, d savoir la determination \ 

des parametres \^t q pour realiser une con^ion d'obliquite NMO ; 

- la figure 8b est un organigramme representant les etapes d'un second 
mode de realisation particulier de I'invention, d savoir la determination 
des parametres V et /? pour realiser une migration PSTM ; 

25 - la figure 9 est uri schema illustrant les difl^6rentes operations realis6es 
afin de detemniner les parametres V et /? permettant d'effectuer une 
migration PSTM. 

Le procede selon I'invention est de maniere generaie un precede de 
traitement des enregistrements de traces sismlques k deport variable dans 
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lequel, d partir desdites traces sismiques enregistr6es, on constitue des 
collections de traces en point milieu commun (CMP) et on soumet las 
traces de chacune des collections a une cpn-ection d'obliquite. 

Le proc6d§ selon I'lnvention detennine en particuller les parametres 
5 de Vitesse V et d'anellipticite n permettant de realiser un traitement 
comprenant une telle correction d'obliquite des traces sismiques d'une 
collection CMP. 

Ledit traitement peut par exemple §tre : 

• une correction d'obliquite NMO anelliptique des traces sismiques ; 

10 . une migration PSTM anelliptique qui, comme on I'a vu pr6c§demment, 
met conjointement en oeuvre les operations de TZO et de migration 
(on parlera cl-aprds de con-ection d'obliquite PSTM). 
La description ci-apr6s concerne plus sp6cifiquement fa con-ection 
d'obliquite NMO. On comprendra cependant. notamment au regard des 
15 deux modes de realisation particuliers de I'lnvention qui seront detallles par 
la suite, que cette description s'applique egalement a tout traitement 
comprenant une con-ection d'obliquite, et notamment un traitement 
comprenant une correction d'obliquite PSTM. 

Afin de realiser precisement la detemilnation de V et de rj, deux 
20 nouveaux parametres dtn, sont consider6s : 

• To qui represente le temps de propagation S deport nul en 
« coonjonnees hj^erboliques » (cf. figure 1), 



• dtn qui represente la correction d'obliquite pour le d6port Xmax le plus 
important (cf. figure 1), 




1+877 



Equation (2) 



dtn = t. 



soit : 




Equation (3) 
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il est important de noter que dtn est dSfini relativement d la vitesse V 
et a ranellipticit^ rj , alors que Tq est un paramStre parfaitement anellipfique 
d§fini relativement d ranellipticit^ 7 , Ind^pendant de V. 

De la sorte, les parametres de vitesse V et d'anellipticite 7 peuvent 
etre calcules, conformement aux equations (2) et (3), en utillsant les 
equations de conversion (4) et (5) suivantes : 

V = Equation (4) 



et, 7j = Equation (5) 

La figure 1a repr^sente Fhyperbote anelliptique d6cal6e permettant 
10 de rSaliser la correction d'obliquite NMO et iflustre la signification des 
parametres et dtn. 

Lesdites « coordonn6es hyperboliques » sont representees sur cette 
figure 1. Leur origine sur I'axe des temps est prise a Tintersection dudit axe 
des temps avec I'asymptote tangente k ladite hyperbole d^calee (cf. 
15 equation (1 a)) a large offset 

En utilisant les deux parametres (dtn, r^), I'equation (1a) de 
riiyperbole anelliptique decaiee devient : 

t = to - CT^^^^:^ Equation (6a) 

V ■^max 

Les parametres (dtn. Tq) definis relativement ^ la Vitesse V S 

20 I'anelllpticite rj permettent ainsi de rendre la conrection d'obliquite NMO. 

CORRnmo = t - to devant dtre appliqu6e aiix traces de deport x 
independante de to : 

CORRnmo(x) = -■ ro + Equation (7a) 

II s'agit done d'une correction d'obliquite statique. Autrement dit, les 
25 donnees enregistrees sur une trace de deport donne seront toutes 
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t = ta-T„+. 



corrlg^es de la mSme mani6re pour un couple (dtn. independamment 
du temps auquel ces donn§es ont 6te acqulses. 

Par consequent. I'estlmatlon des param^tres de Vitesse et 
d'anellipticit6 n'est pas perturb6e par r6tirement (« stretch » selon la 
terminologie anglo-saxonne) des traces generalement observe lorsque des 
con-ections d'obIlquit§ dynamlques sent r6alisees. 

De mani^re similaire. en utillsant les deux parametres (dtn.r„). 
r^quation (1 b) de la double hyperbole d^calee DSQR de la migration PSTM 
aneiliptique devient : 

.MI^^i^SIE^^ 6 . dta(dtp.^2.„yx3Z: :^ ^ 

V 4 ^VT"-^ ^ ^ Equation (6b) 

La figure 1b represents la double hyperbole d6cal6e DSQR de la 
migration PSTM parametree avec les paramdtres (dtn. rj. 

Dans le cadre de la migration PSTM. x^, represents le maximum de 
deport et d'ouverture de la migration. 

Les parametres (dtn.r^) d6finis relativement a la vitesse V et ^ 
ranellipticite rj pemiettent ainsi de rendre la correction d'obliquite PSTM 
CORRpsTM = t - to devant etre appliquee aux traces de deport x 
independante de to (equation (7b)) : 

II s'agit done ^El'une con^fon PSTM statique. AutremTnt dit 
rensemble des echantillons d'ouverture x-x„, d'un cube de traces ^ deporl 
constant (cube « iso^ffset ») est. pour un couple (dtn.rj donne. decaie 
d'un m§me temps. 

En analysant les parametres (dtn.rj en une pluralite de temps de 
pointe. le procede selon I'invention permet en particulier de determiner les 
parametres (V,,;) necessaires pour realiser un traitement c;)mprenant une 
con-ection d'obliquite aneiliptique des traces d'une collection CMP 
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Ledit precede comporte de manidre simpiifi^e les 6tapes pr6sent§es 
ct-aprSs. 

Au cours d'une ^tape pr^liminaire, un volume d'anaiyse comprenant 
une plurality de nceuds (dtn, Tq) est dSfini. 
5 Puis, pour chacun des noeuds de ce volume, on effectue : 

• dans un premier temps, selon l'6quation (7a), la correction d'obliquit^ 
statique des traces de la collection CMP etudiee, en fonction des valours 
des parametres dtn, au noeud consid^r^, ladite correction statique 
6tant valable pour n'importe quel temps du pofnte ; 

10 • dans un second temps, on calcuie fa semblance, fonction du temps, 
associSe a la correction r§alis§e i I'etape pr§c§dente. 

Enfin, pour ciiaque temps de la pluralite de temps de points, on 
determine : 

• dans un premier temps, le noeud (dtn, t^) permettant une correction 
15 optimale, par exemple au regard du crit^re de semblance (un tel crit^re 

§tant ciassiquement utilise en traitement sismique pour « mesurer 
rhorizontalite » des courbes de reflexion et determiner la fiabilite du 
pointe) ; 

• dans un second temps, on convertit les valeurs des parametres dtn, 

20 audit noeud de semblance maximale en valeurs des paramdtres de 
Vitesse V et d'anellipticite rj audit temps de pointe consider^. 

La lol de vites^ (c'est-a-dire I'ensemble des couples (temps du 
pointe, V)) et la lol d'anellipticite (ensemble des couples (temps du pointe, 
n)) sont ainsi etablies. 
25 Finalement la correction d'obliquite de I'ensemble des traces 

sismiques peut etre realisee en utilisant ces lois de Vitesse V et 
d'anellipticite n dans I'equation (1a) de I'hyperbole anelliptlque decalee. 

De maniere similaire, et comme cela sera plus particulierement 
detailie par la suite, les parametres (dtn.ro) 'a migratiop PSTM statique 
30 (cf. equation (7b)) peuvent egalement §tre pointes. Les lois de vitesse V et 
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d'anelliptrcite n sont alons d6termln6es et peuvent gtre utllis§es. pour 
r§aliser la migration PSTM. dans I'equation (1b) DSQR de la double 
hyperbole anelliptique decalee. 

La figure 2 lllustre I'effet du paramdtre r„ (et done. cf. Equation (2). 
de ranellipticlte tj) sur ies r6sldus de courbure apr^ correction d'obliqult^ 
NMO. 

On notera que I'^chelle verticale tennporelle de la courbe de la figure 
2 est exager6e afin de bien rendre compte de cet effet. 

Trols courbes sont representees sur la figure 2 pour lesquelles le 
param§tre dtn est fix§ ^ la valeur correcte et le paramdtre r, prend 
diff§rentes valeurs. 

La courbe centrale repr6sente le cas oCi est ^ sa valeur correcte 
To,, c'est-^-dlre lorsque I'anellipticite 7, con-espondante est S sa valeur 
con-ecte 7^,. Comma cela est attendu, la courbe de reflexion comgee est 
alors horizontale. 

La courbe sup6rieure repr6sente le cas oO t, prend une valeur 
inf^rleure d sa valeur correcte ^r„.. I'anelliptlcitd rj, correspondante etant 
sup§rieure d sa valeur con-ecte tj 

I true ' 

La courbe inf6rieure repr6sente quant ^ eile le cas oO prend une 
valeur v,, sup6rieure d sa valeur correcte r,^, ranellipticlte 7, 6tant ' 
inferieure a sa valeur correcte n 

/ true ' 

On remarque de ces courbes inferieure et sup6rleure que 
l'« horlzontalite » de la courbe de reflexion con-|g6e est acceptable a faible 
deport (X «0) et a large deport (x «Xn,ax). 

En revanche, on observe des residus de courbure signiflcatlfs 
lorsque le deport x ne tend pas vers une de ces valeurs limites 0 et x„ax. On 
note en particulier des r6sidus particulierement signlficatifs pour un deport x 
centre au milieu de la gamme de deports. 
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A titre d'exemple, et comme cela apparalt sur la figure 2, lorsque 
est li ro2, une correction r6siduelle RMOto^ devrait %tre apportee aux 
traces de deport x^^^,^ . De manidre simiiaire, lorsque est a r^j, une 

conrectlon rSsiduelle RMOtoj devrait §tre apportee aux traces de deport 

Du fait de ces r^sidus de coiirisure significaOfs. la gamme des 
deports peut quasiment Stre utilis§e dans son int^gralite afin de determiner 
rapellipticite ij. 

Le param^trage en (dtn, r^) de la correction d'obliquite rend done 
10 possible rutilisation des donnees disponibles pour I'ensemble des deports 
(x compris entre 0 et Xmax) dans la determination de ranellipticlte tj. 

Or, comme cela a dejS 6t§ evoqu^ pr^cedemment, cela n'est pas le 
cas des corrections d'obliquite paramStrees avec V et 77 pour iesquelles 
sont essentlellement les donnSes a large d§port qui permettent restimatlon ' v 
15 de ranellipticite 77 / 

Ainsi, Teffet du nouveau parametre anelliptique Tq est distribue sur v 
tous les deports, au contraire de ranellipticit6 tj qui n'affecte que les larges 
deports. Le « comportement » du paramdtre permet done de mieux 
contraindre les valeurs d'anelllpticit^. . \^ 

20 Comme cela a 4r§ja et6 mentlonn§ pr6cedemment, la determination 

des parametres dtn et optimum est r^allsee dans un volume d'analyse 

3D(to, dtn. To). 

J* 

On considere dans ledit volume d'analyse une plurality de noeuds 
(dtn, To), c*est-a-dlre une pluralite de couples de parametres dtn, dont 
25 les valeurs sont connues. 

Les noeuds sont g6n6ralement espac6s r§gulierement les uns des 
autres, d'un increment Adtn sur I'axe dtn et d'un increment' Aro sur Taxe . 



1erd6p6t 



16 



10 



Des valeurs minimales dtn^,„. t,^,„, ton„n et maximales dtn^. 
tomax des parametres respectlvement dtn. et to permettent de d^finfr les 
limites dudit volume d'anaiyse. 

De manidre avantageuse, des valeurs plausrbles des paramdtres de 
Vitesse V et d'anellipticit§ 17 peuvent §tre utillsees afin de d§finir S I'int^rieur 
dudit volume d'anaiyse un corridor [dtn.,n(to). dtn.3x(to)].[r„.i„ (to), r,.3x(to)]. 

Ce corridor permet de restrelndre le volume d'anaiyse et done le 
nombre de noeuds (dtn. r,) devant §tre consid6r6 pour la determination du 
couple (dtn, To) optimal. 

SI I'utilisation dudit corridor est b6n§fique pour I'efficacite du precede 
selon I'inventlon. elle pemiet §galement de contraindre la solution vers les 
bons ph6nom§nes, sans avoir ^ considerer des couples (dtn, (et done 
des couples (V, 77)) incompatibles. relatifs par exemple k des reflexions 
multiples ou & divers ph6nom6nes d'interfference. 

Les figures 3a et 3b repr6sentent ledit volume d'anaiyse (to, dtn, rj 
dans le cadre d'un exemple de traitement de donnees sismiques reelles 
lors d'une connection d'obliquit6 NMO r^alisee selon le precede de 
rinvention. 

La figure 3a repr§sente trois panneaux 2D a, b et c du volume 
20 d'anaiyse : 

- le panneau ^ est un panneau (dtn, S to constant ; 

- le panneau b est un panneau (dtn, to ) S t,, constant ; 

- le panneau c est un panneau (^o . to ) a dtn constant. 

La figure 3b repr§sente de maniere sch^matique le volume d'anaiyse 
25 3D (to. dtn, rj ainsi que trois intersections a ce volume selon trors plans a 
respectivement to. t„ dtn constant, chacune de ces intersections 6tant 
projet6e sur le panneau a, b ou c correspondant de la figure 3a. 



15 



17 



Les panneaux (dtn, ro) a to constant (panneau a dans l*exemple ci- 
dessus) sont les panneaux dans lesquels est r6alis6 le pointe bispectral des 
paramdtres dtn et r^, par exemple selon le critdre de semblance maximale 
de la correction d'obliquit^, pour le temps de point6 to considere. 

Le corridor [dtnmin(to), dtnmax(to)] I^omin (to), ^omaxCto)] mentlonne 
pr6c6demment pour {'analyse effective k f Int6rieur du volume d'analyse est 
§galement repr6sente sur ia figure 3a (zone plus claire sur chacun des 
panneaux). 

Le pas d*6chantillonnage en temps Ar^ des donn§es sismiques 
definit I'ecart entre deux temps de points successifs (temps pour lesquels 
on determine le nceud (dtn, Tq) de semblance maximale associ§) et done le 
nombre de panneaux de pointe bispectral devant §tre consid§r6. 

Le pointe automatique permet ainsi d'extraire ies param^tres dtn, 

en une denstte d'autant plus importante que Tincrement Ato entre les temps 
de points est faible. 

L'echantillonnage des parametres d'anaiyse dtn et est quant a lui 

directement lie a la resolution de la sismique, dtn et ayant effectivement 
les m§mes dimensions que les enregistrements sismiques. 

La recherche systematique du maximum de semblance, 
classiquement connue en soi, permet de determiner la paire (dtn, r^) 
procurant, pour un terras de pointe to donne, la meilleure focalisation. 

Des interpolations paraboliques autour des valeurs des noeuds (dtn, 
To) peuvent en outre permettre d'evaluer les valeurs des parametres dtn, 
To entre les differents noeuds effectivement point6s. Et une telle evaluation 
rend en particulier possible la d6termination plus precise encore (par 
opposition §1 la determination limitee aux noeuds du corridor) du couple de 
parametres dtn, maximisant la fonction de semblance. 
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Flnalement les param6tres de Vitesse V et d'anelllpticit6 rj sont 
d6termin6s, toujours pour le temps du points to consld6re, en utilisant les 
Equations de conversion (4) et (5) susmentionnees. 

La figure 4 permet de comparer les deux approches (V, Van) et (dtn, 
rJ en repr6sentant, pour un temps de polnt§ donn6, leur panneau de 
point§ bispectiBi respectif. 

La figure de droite illustre approche classlque (V. Van) pour laquelle 
les deux axes sont des axes de Vitesse (ranellipticlte rj §tant lie au rapport 

de ces deux vitesses selon ij = - 



ry 4 N 
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La figure de gauche illustre quant d elle I'approche (dtn. x^) selon 
I'inventlon pour laquelle les deux axes sont des axes temporels. 

II est Important de noter de I'etude de cette figure 4 que les 
parametres dtn et r^, semblent etre d6corr6l§s. Cette « d6corr6lation » est 
frappante lorsque I'on compare les deux approches. I'^talement du spectre 
(dtn. To) etant effectlvement beaucoup plus r6durt que celui du spectre (V, 

Van). 

Le points r§alls6 dans le cadre de I'approche (dtn. r^) selon 
{'Invention est done plus precis que celui classiquement realise. 

En outre, cette « d^correlation » permet le filtrage des point^s dtn et \. 
separement, tout ^ conservant les corrections d'obliquit6. Cela n'est 
pas le cas pour les parametres V et /; pour lesquels la diminution de I'un 
des parametres doit imp§rativement §tre compensee par I'augmentation de 
I'autre, et vice versa. 

Et grSce aux Interpolations et filtrages individuels des parametres de 
rinvention. dtn et ^o. il est alors possible de r^aliser I'interpolation et le 

filtrage simultan^s des paramdtres standards V et /? de correction 
d'obliquite. 
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La figure 5 reprSsente la correspondance non linSaire, selon les 
Equations (4) et (5) presentees ci-dessus, entre la paire de param^tres 
temporels (dtn, r^) et ia paire de param§tre de vitesse (V, Van). 
La figure 6 represente de gauche a droite : 
5 - une collection CMP de traces sismlques reelles avant correction 
d'obllquite ; 

- le pointe du parametre dtn correspondant k cette collection de traces ; 

- le pointe du parametre correspondant a cette collection de traces. 

Sur le pointe du parametre § droite sur la figure 6, la ligne droite 
10 ro = to correspond aux courbes de reflexion purement hyperboliques. 

La figure 7 represente quant a elle, de droite S gauche, la fonction de 
semblance et les fonctions de vitesse V et d'anellipticite rj deduites (cf. 
equations (4) et (6)) du point§ des paramdtres dtn et represents sur la 
figure 6. 

15 On note, sur cet exemple de traltement de donnSes sismlques 

reelles, que les valeurs de V, rj obtenues correspondent pour la plupart a 
une semblance superieure a 40%. 

La description ci-apres detaille deux modes de realisation partlculiers 
de rinvention. 

20 Le premier de ces modes concerne un precede de determination des 

parametres optimum pour realiser une correction d'obliquit§ NMO 

.it 

anelliptique des traces d'une collection CMP (voir les differentes etapes 
representees sur rorganigramme de la figure 8a). 

En reference S la figure 8a. ce premier mode de ^ realisation 
25 comprend une etape 1a d'Initlalisation au cours de laquelle sont realisees 
successlvement les operations de : 

• determination des limites du volume d'analyse [dtnmin, dit^maMu 

mini ^0 maxlfPomfnt tomaj I 

f 
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. calcul des corrections d'obliquite CORR^mo (equation (7a)) pour 
tous les deports et pour tous les noeuds (dtn. inclus dans le 
volunfie d'analyse ; 

• delimitation ^ Flnterieur du volume d'analyse du corridor [dtn^„(t6). 
5 cItnmax(to)] [ronf,in(tb), T^tMio)] des valeurs plausibles de vltesse et 

d'aneiliptlcit§ ; 

Une fols cette 6tape la d'inltialisation r^ailsee, une §tape 2a de 
calcul des lois de Vitesse V(to) et d'anelliptlclte n(to) est mise en oeuvre pour 
chaque collection de traces CMP. 
10 Cette etape 2a comprend : 

. une premiere operation 3a r^alisee pour chaque noeud (dtn. r,) du 
corridor d^fini lore de r§tape la d'inltialisation. au coure de laquelle 
sont realisees successivement. pour chaque temps du pointe to, les 
operations de : 

- application pour tous les deport, le long du corridor, des 
connections CORRnmo d'obIiquit§ statiques pre-calcul§es lors 
de l'§tape la d'initialisation ; 

- calcul de la fonctlon de semblance sur les donn6es corrigees 
le long du conidor en utillsant une fenetre temporelle 

^ appropriee a i'ondelette dominante ; 

- sommatlon (calcul du « stack ») des donnees corrigees le ' 
long du^corridor (seules les donnees a faibles d§ports 
pouvant avantageusement §tre pour cela utilisees) ; 

. une seconde operation 4a r6alis6e pour chaque temps tb'du polnt6 
(lesdits temps itant espac6s de Ato entre [to„un. ic^l au coure de 
laquelle sont r6alis6es les operations de : 

- recherche du maximum de semblance dans le corridor 
[dtn„,in(to). dtn„^(to)],[Tomi„(to). T,„,^(to)] du panneau 
bispectral (dtn, r^); 
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- cqntrdle du fait que la position en (dtn, r,,) du maximum de 
semblance conrespond d un extremum de la sommation pour 
les m§mes vaieurs dtn et To ; 

- creation des series dtn(to), (to) et semblance (to) ; 

• une troisiSme operation 5a visant a sSlectionner et ajuster les point§s 
obtenus, au cours de laquelle sont r§alis6es les operations de : 

- tri croissant de la s6rie semblance(to) ; 

- validation des pointes dtn et pour lesquels la distance en 

temps aux pointes de semblance plus 6leve est superieure a 
une valeur pr§d6finie ; 

- ajustement des vaieurs dtn et polnt6es et validees par des 
interpolations paraboliques utilisant des vaieurs alentours ; 

- retention des vaieurs pointSes, valid6es et ajust6es si le 
calcul des vitesses d'intervalle de Dix avec les point6s de 
semblance plus Slev^e est possible. 

• une quatrieme operation 6a visant S convertir, ^ I'aide des Equations 
(3) et (4), les vaieurs de dtn et pointees, validees, ajustees et 

retenues au cours de I'operation 5a en lois de Vitesse V et 
d'anellipticite n- 

Les lois fonctions du temps de la Vitesse V et de ranellipticit^ rj sont 
ainsi parfaitement det^in6es. Et la correction d'obliquit6 NMO anelliptique 
des traces sismlques de la collection CMP peut ainsi dtre pr§cis6ment 
MalisSe. 

Le second mode de realisation particulier de I'lnvention ^jnceme un 
precede de determination des paramStres optimum pour la migration PSTM 
anelliptique des traces d'une collection CMP. 

Ce second mode de. realisation peut §tre comprls comme une 
generalisation du premier mode discute ci-dessus. 
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En effet. comme cela a 6te montre pr6c6demment. rutllisatlon des 
parametres (dtn. rj permet des corrections d'obliqult§ PSTM statiques (cf. 
Equation (7b)). 

Cette utilisation pr§sente. dans le cadre de la con-ectlon d'obliqurte 
PSTM, les m§mes avantages que ceux pi^cedemment discut6s pour la 
correction d'obliquite NMO. 

Plus preclsement, on notera que le premier mode de realisation n'est 
qu'un cas particuller du second mode de realisation con-espondant au cas 
d'une ouverture de migration nulle. 

En reference ^ la figure 8b. ledit second mode de realisation 
comprend une 6tape 1b d'initialisatlon au cours de laquelle sent realis^es 
successivement les operations de : 

. detemilnation des limites du volume d'analyse[dtn„.^. dtn^.fr, 

^0 maJ,[tomfnf tomaxl I 

• calcul des conections d'obllqult6 CORRpstm (equation (7b)) pour 
tous les noeuds (dtn, r„) indus dans le volume d'analyse et pour 

tous les deports de migration S rint6rieur de Touverture de 

migration ; 

• delimitation ^ I'int^rieur du volume d'analyse du conidor [dtn^h(fo). 
citnmax(to)],[ron,in(to), z-omax(t6)] des valeurs plausibles de vitesses et 
d'anellipticite ; • V 

Une fols cette 6fepe 1b d'initialisatlon r6alis6e. ledit premier mode de 
realisation met en oeuvre. pour chaque collection de traces CMP, une etape 
2b de calcul des lois de vitesse V(to) et d'anellipticite n{to). 
Cette 6tape 2b comprend : 
• une premiere operation 3b r6alisee pour chaque noeud (dtn, r„) du 
corridor defini lors de I'etape 1a d'initlalisation. au cours de laquelle 
sent realis6es successivement : 

o pour chaque classe de deport, les operations de : 
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" application sur tous les points milieu a I'interieur de 
Touverture de la migration, le long du comdor, s 
confections CORRpstm statiques prS-caiculees lots de 
rstapelbd'lnitialisation; 

■ sommation des points milieu corrigSs le long du 
comdor ; 

o pour chaque temps du points to, les operations de : 

■ calcul de la fonction de semblance sur les donnees 
corrig6es le long du corridor en utilisant une fenetre 
temporelle appropriee a Tondeletle dominante ; 

■ sommation (calcul du « stack ») des donnees 
corrig^es le long du corridor (seules les donnees d 
faibles deports pouvant avantageusement etre pour 
cela utilisees) ; 

• une seconde operation 4b r^alisde pour chaque temps to du points, 
similaire a I'operation 4a decrite pr^cedemment, visant a cr^er les 
series dtn(to), ro(to) et semblance (to) ; 

• une trolsi^me operation 5b, similaire a reparation 5a decrite 
precedemment,visant a s61ectionner et ajuster les pointes obtenus. 

• une quatrieme operation 6b, similaire a I'operation 5a decrite 
pr6c6demment,visant ^ cohvertlr les valours de dtn et Tq en lois de 
Vitesse V et d'an^ipticit6 n- 

Les lois fonctions du temps de la vitesse V et de Tanellipticitd rj sont 
ainsi parfaitement d6termin6es. Et la migration PSTM des traces sismiques 
de la collection CMP peut ainsi etre precisement r^alisee. 

La figure 9 illustre le second mode de realisation de Tinventfon dont il 
vient d'etre fait 6tat 

Les donnees sismiques sont initialement regroup^es en cubes a 
deport constant (cube iso-deport ou « iso-offset »). 

Comme cela vient d'etre presente, les corrections d'obliquite 
CORRpstm sont appliquees (cf. premiere operation de Tetape 3b), pour 
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chaque classe de d6port (soit pour chaque cube iso-deport). sur tous les 
points milieu. La fleche llbellee MIG sur la figure 9 illustre cette operation. 

Les points milieux ainsi corriges sont ensuite somm6s au cours de la 
seconde operation de I'etape 3b. la fleche Iibell6e STACK sur la figure 9 
illustrant cette operation. 

Ces deux operations MIG et STACK sont sp6cifiques au second 
mod© de realisation de Tinventlon (migration PSTM). 

Les operations sulvantes sont en revanche mlses en oeuvre dans le 
cadre de chacun des deux modes de realisation discutes (la con-ectlon 
d'obliqulte NMO correspondant. comme on I'a vu. au cas particulier de la 
migration PSTM d'ouverture de migration nulle). 

Pour chaque temps du pointe. la semblance est ensuite calcuiee 
(seconde operation de I'etape 3a pour la seule correction d'obliquite 
troisidme operation de I'etape 3b pour la migration PSTM). la fleche libellee 
•semblance' sur la figure 9 illustrant cette operation. 

Bien entendu, les operation MIG. STACK et de calcul de la 
semblance sont mises en oeuvre pour chaque nceud (dtn.r^). 

Le 'pointe bispectral automatique' mentionne sur la figure 9 
con-espond au pointe des parametres (dtn.r„)de semblance maximale pour 
chaque temps todu pointe (operations 4a. 4b sur les figures respectivement 
8a. 8b). La fleche llbeliee MAX illustre la recherche du maximum de > 
semblance. ^ 

Finalement les pointes des parametres (dtn.rj pointes sont 
convertis en lois de vitesse V et d'anellipticite (operation 6a. 6b sur les 
figures 8a, 8b et fleche libeliee CONV sur la figure 9). 



25 



REVENDICATIONS 

1. Proc6d6 de determination des paramfetres de vltesse V et 
d'aneliipticite n pour un traitement de traces sismiques d'une collection 
h point-milieu commun {CUP) comprenant une correction d'obliquite 
anelliptique, caracteris6 en ce qu'il comporte : 

t une etape preliminaire de definition d'une plurality de noeuds 
(dtn, ), lesdits noeuds §tant significatifs de paramdtres dtn et 
repr^sentant respectivement ia correction d'obliqult6 pour le d6port 
maximal et le temps de propagation S deport nul en coordonn^es 
hyperboliques, iadite §tape preliminaire etant suivie 

• pour cliacun des noeuds (dtn, r^) d^finis lors de I'Stape 
preliminaire, des etapes de : 

- correction d'obliquite statique des traces de la collection 
CMP en fonction des valeurs desdits parametres dtn et 
au noeud considere, et de 

- caicul de la fonction de semblance associee a Iadite 
correction d'obliquite p>our le noeud considere : et 

• pour chaque temps to du pointe, d'une etape comprenant la 
determination du noeud (dtn(to),ro (to)) de semblance maximale, 

• et d'une etipe finale de conversion des parametres dtn(to) et 
ro(to) determines pour chaque temps to du pointe de maniere k 
obtenir les lois de Vitesse V(to) et d'anellipticite /7(to). 



ler d6p6t 
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2. Proc6d6 selon la revendlcation 1 . caract§rls6 en ce que les noeuds 
sont d§finis lors de I'etape preliminaire dans un volume d'analyse 
(dtn,iro.to) determine par des valeurs minlmales et maximales 

respectlvement [dtn,rtn.dtn,aJ.[To„„n.romax] et [tdn,.„,t6„,aj des paramdtres 
dtn, Tf, etto. 

3. Proc§d6 selon la revendlcation 2, caract6rise en ce que, au cours 
^® i'etape pr6limlnaire, un con-ldor 
ldtnmm(to),dtnmax(to)].[romin{to).roniax(to)] d*6volutlon des param^tres dtn et 
To est delimits d I'lnt6rieur du volume d'analyse en fonction de valeurs 
plausibles des param6tres de vftesse V et d'anellipticite rj, les nceuds 
(dtn. z-o) d§finis pour rappHcation de la correction d'obliqulte etant alors 
situ6s le long du comdor ainsi d61imlt§. 

4. Precede selon Tune des revendicatlons prec6dentes. caract6ris6 en 
ce que. pour chacun des noeuds (dtn.r^). 11 comporte en outre, suite a 
r§tape de calcul de la fonction de semblance, une etape de sommation 
des traces sismiques corrig^es. 

5. Proc6d6 selon la revendlcation 4, caract§ris6 en ce que la 
sommation des t^^ces comgees est realisee en n'utilisant que les 
traces d faibles deports. 

6. Proc6d6 selon I'une des revendlcations 4 ou 5, caract§ris6 en ce 
que. pour chaque temps de points, il comporte en outre. suite ^ I'dtape 
de determination du noeud de semblance maximale, une etape visant i 
controler que les valeurs dtn et du noeud de semblance maximale 
con-espondent ^ un extremum de la sommation pour les memes valeurs 



dtn et To . 
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7. Proc6d§ selon I'une des revendications pr6cedentes. caracteris6 en 
ce qu*H comporte en outre, suite a I'etape mise en oeuvre pour chaque 
temps to de pointe pour la determination du noeud (dtn{to), r^ito)) de 

semblance maximale, et prealablement a Tetape de conversion, une 
etape de selection et d'ajustement des pointes effectues. 

8. Proc§d6 selon la revendlcation 7, caract6ris§ en ce que ladite 6tape 
de selection et d'ajustement des pointes comporte une etape visant ^ 
ne conserver que les point§s dtn et pour lesquels la distance en 

temps aux point6s de semblance plus Sieve est superieure d une valeur 
pr6definie. 

9. Procede selon la revendication 8, caract6rise en ce que ladite etape 
de selection et d'ajustement des pointes comporte en outre une etape 
visant a ajuster les pointes conserves dtn et par interpolations 
paraboliques utilisant des valeurs autour desdites valours pointees. 

10. Precede selon la revendication 9, caractSrise en ce que ladite Stape 
de selection et d'ajustement des pointes comporte en outre une §tape 
visant a eliminer les pointes dtn et conserves et ajustSs lorsque le 

calcul des vitesses d'intervalle de Dix entre le pointe considere et les 
pointes de semblance plus elevee est impossible. 

11. Proc6d6 selon Tune des revendications 1 & 10, caract6ris6 en ce 
que le tra'rtement applique aux traces sismiques est un traitement de 
correction d'obliqulte NMO mettant en oeuvre une correction d'obliquite 
statique CORRnmo. 
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12. Proc6d6 selon la revendicatlon 11 , caract6ris§ en ce que. au cours 
de l'§tape prelimrnaire. les corrections d'obliquite CORRnmo sont 
calculees pour tous les nceuds (dtn. r,) Indus dans le volume d'analyse 
et tous les deports des traces slsmlques tralt6es. 

13. Proc6d6 selon la revendicatlon 12. caract6ris6 en ce que I'^tape de 
correction d'obllqult§ r6alls6e pour chacun des noeuds (dtn. r^) 

conslste S appllquer les corrections d'obliquite CORRnmo calculees lors 
de l'6tape pr6limlnaire. 

14. Precede selon I'une des revendications 11 ^ 13. caract§ris6 en ce 
que, pour un couple (dtn, t,) donn§. la correction d'obliquite statique 

CORRnmo d'une trace sismlque de d§port x est r6alis6e selon 
I'equation : 

repr6sente le deport maximal de la collection CMP. 

15. Precede selon I'une des revendications 1 i 10. caract6ris6 en ce 
que le traitement appllqud aux traces sismiques est une migration 
PSTM mettant en oeuvre une con^ctfon d'obliqult6 statique CORRrstm- 

16. Proc6d6 selon la revendicatlon 15. caract6rls6 en ce que, au cours 
de l'6tape prelimlnaire. les con-ections d'obliquite CORRpsrw sont 
calculees pour tous les noeuds (dtn. r,) Indus dans le volume d'analyse 

et tous les deports de migration ^ I'lnt^rieur de I'ouverture de la 
migration. 
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17. Proc§d§ selon la revendication 16, caracteris^ en ce que I'etape de 
correction d'obliquite realisee pour chacun des nceuds (dtn, r^) 
comporte, pour chaque classe de deport, I'application des corrections 
d'obliquite CORRrstm calcul6es lors de r§tape preliminaire sur tous les 
points milieu a Tinterieur de Touverture de la migration. 

18. Procede selon la revendication 17, caracterise en ce que I'etape de 
correction d'obliquite realis6e pour chacun des nceuds (dtn, t^) 

comporte en outre, pour chaque classe de deport, suite d ladite 
application des conrecHons d'obliquite CORIRrstm, la sommation des 
points milieux corriges, 

19. Proc6d6 selon Tune des revendication 15 gi 18, caract§ris6 en ce 
que, pour un couple (dtn, r^) donn§, la correction d'obliquit6 statique • 
CORRpsTwest realisee selon Tequatlon : 




dans iaquelle : 

- Xm represente les coordonnees des points milieux, 

- x-Xm repr6sente Touverture de la migration PSTM, 

- h repr§s^te le demi-d6port source — recepteur, 

- Xmax represente le maximum de dSport et d'ouverture de la 
migration. 

20. Proc6d6 selon Tune des revendications 14 ou 19, caract§ris6 en ce 
que, au cours de l'6tape finale de conversion, les param^tres dtn(tb) et 
To (to) sont convertis en loi de vitesse V(to) selon T^quation 
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21. Proced§ selon I'une des revendications 14 ou 19, caracteris6 en ce 
que, cours de l'§tape finale de conversion, le param^tre To(to) est 

converti en loi d'anelliptlcit6 nih) selon n = -( -^-i^ 

22. Proc6d6 selon les revendications 20 et 21, caract6ris6 en ce que le 
param^tre dtn est d^fini relativement § la Vitesse V et ^ I'anellipticite rj 
selon r^uation : 

dtn = --^c, + 

23. Proc6d6 selon la revendlcation 21, caract6ris6 en ce que le 
paramdtre t„ est d6fini relativement S ranelliptfclte rj selon I'equation 

" 1 + 8^- 

24. Proc§de pour caracteriser un champ de vitesses en vue d'un 
traitement de donnees sismiques a partir d'une collection ^ polnt-mflieu 
commun (CMP) de traces sismiques. caracteris§ en ce que I'on definit, 
pour chaque temps to de propagation pour un d§port nul. un ensemble 
de param^tres dtn et repr§sentant respectivement la con-ecfion 
d*obliquit6 pour le deport maximal et le temps de propagation a deport 
nul en coordonn^es hyperboliques. 
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D<:lcci«ijaatioudes limitejs du volurab d'cmalyse (dtii> to, to) 
Calcul des coireclions d'oblk|intc CORRkkto pom* tou^ les 
d6poit« et pour tous tes ngeuds (din, to) 
Calcul.d«s liraitw du corridor des analyjjes effectives 




Pour clj^quc le loug d« liwrf drtr: 

^ Appitcatton descoirections d'obHqiut6 CORRkmo 

• QhIoiiI di? la sciublanq^ 

• Calcui dc la sominati«iii eti d6ponspro*6hes 



l^cctctdiQ du m?»?cimx«if) 4^ sem'blarice dkns le corritlor 
ct iicend.{dtii^ To)-'cc)>fT«§p<5hdant 




'Tri Groissanticle la- seri'e sertih1arice(fe) 

R^gej d.e5 poiijtte& JUtjp. p^^p5&iics dfes^^ -j^jQijnt^s i plus forte' 

Ajustement par int«jpolatloiis paraboliques d^ valeurs 

Tfeejet des couples (din, To).doiit Ics vitesses d^otemlle de 
JM ayeclesWuplgs (dtn, t(>)'dpjit la seniblanc^^ est i)lus 
forle^ soni inacccpvablcs. 
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CalculdeshmitesducomdorL^U effectives 



Pow chaque t« le long da confdor: 

Apphcation des corrections d'obUauite ropu 
• Calculde la semblance ° CORR^mo 
Calcul de la sommation en deports proches 



et du ncBud (dto, Tt>) coirespon dant 



Lenioeud(dtn,To) 



conespond-dl4un 
extremum du stack ? 



Tri WMssant de la serie semblance(to) 
fSilto.'^'^ PO-t^ ^ Pl- forte 

^ir;TSe4:,^^^'^°^^^^ ~<^- valeu. 

foite, so^tt ^ -"Wance est plus 





Deterniinalion dcs liniilcs du volume d* analyse (dta, to, lo) 

• Calwl de5 correcrit^jis d'obliquile CORRpi^m pour Imiii les 
nocuds (dtii, To) fiss deports & rinlerieur do ronve:iTi?re de 
migration 

• Calcul dbs limites da corridor "des analyses cffeciivcs 




POU r chaque ctal;se de deport 

► AppUeation sur tons les points milieux a riatdneur de 
rouverture de la migralion, Iclong du conidor, des com^ciiuns 
d'obIj.qiiiteCOkRp<;T^i 

• S.pmmalioii dp^ points milicmx conrig^ Ve long du corridor 
Pour cjiaque to le long du carridjur: 

Galcid du la semblance 

* Calcul db la isommatioii §n depoHs proches 




Recherche- du ms^cimirm dc semblaitce djrnis le ooiridpr 
el dit'xpjud (din, to) correspondaiit 






Tri crotssaint de la sede se^iiblanoe(l^i) 
ReJpL.a^.poiVitSs ti-op -^i-ocb^s dis poinxt^ A . "plus. foiTc 

Aju$i94ii^ttt par fnteijix^Iaripns j^ar^tfoliques- des valours 
(din, To) Tolonuos 

Rejeil ^es cipiflp|!csvi(dtn- to) dpnt les \'itessjss d^intervt^lle 
de .Dix avec couples (dto, Toj dotitla s^ata^ 
plus iQric, :sont iffiicc.eptabfes. 
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D^terminatic 




^ourdiaqae de^ipSrt ~^-r^.-^^^^.-Tr...«»««,^i«^^ 

• Calculde la semblance 

Calcul dc Ja sommation en dgporte ryr..i^ 



etdunoeud(dtD,To)coirespondant 
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,ggatio ndes series dta(to), tb(t^ ets^la»cert.^ 



Tri croissant de la serie semblance(to) 
s'tS^r^^*^ ^^^'^^ ^ PI- forte 

Rejet des couples (dtn, to) dont les vitesses d'intervalle 
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